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Abstract 

The annihilation of the Be meson to leptons and inclusive light hadrons is analyzed in the 
framework of nonrelativistic QCD (NRQCD) factorization. We find that the decay mode, which 
CN i escapes from the helicity suppression, contributes a sizable fraction width. According to the 

analysis, the branching ratio due to the contribution from the color-singlet component of the meson 
1^ ■ Be can be of order 10~^. We also estimate the contributions from the color-octet components. 
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With the velocity scaling rule of NRQCD, we find that the color-octet contributions are sizable 



I too, especially, in certain phase space of the annihilation they are greater than (or comparative 

■ to) the color-singlet component. A few observables relevant to the spectrum of charged lepton are 



suggested, that may be used as measurements on the color-octet and color-singlet components in 



Qh' the future Be experiments. A typical long distance contribution in the annihilation is estimated too. 
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I. INTRODUCTION 



The observation of the Be meson by CDF collaboration at Fermilab|l[] is one of the 
important discoveries in heavy quark physics and their results are consistent of the theoretical 
predictions^, ^ ^, ^ 0, when the theoretical uncertainties and experimental errors are 
concerned. More Be events in the new run of the Tevatron than those of the discovery and 
much more events by several orders at the hadron collider LHC are expected^. Therefore, 
thorough experimental study of the Be meson and more precise comparisons with theoretical 
predictions on its properties, especially, on its various decay modes will be available in 
foreseen future. 

One of the interesting decay modes is its pure leptonic decay. Since the decay amplitude 
is proportional to the decay constant fs^ and the CKM matrix element Vcb- In principle, it 
can be used to measure fs^ if we know the value of Vcb, or vice versa. However, the pure 
leptonic decay modes suffer by the the helicity suppression with small leptonic mass. The 
suppression is not severe for Be rvr decay mode only, whereas, this decay mode is hard 
to be detected because its resultant state contains two neutrinos at least. 

There is no helicity suppression if more particles are included in the decays. Such an 
example is the radiative leptonic decay. Be ll'^^u which has been analyzed in p, ^ 
|Tl], 0. The results show that the decay branching ratios are order of 10~^ for / = (/x, e), 
which is much larger than that of the pure leptonic decay of e-channel and comparable to 
that of /i-channel. The relatively larger branching ratio may make them possibly to be 
measured at the LHC. 

In this paper, we analyze a new annihilation decay mode of the meson i?c, i.e., the 'semi- 
inclusive' decay B^ l^vi + light hadrons. Here the light hadrons are produced from the 
emitted soft or/and hard gluons, i.e., we are not interested in any specific light hadron(s) 
but sum all of them. The annihilation process escapes from the the helicity suppression due 
to additional gluons in the final state. Since it is an annihilation of a pair heavy quarks 
essentially, so it can fully be analyzed in the framework of Nonrelativisitic QCD (NRQCD) 
factorization |T^. According to the framework, the decay width can be factored into a sum 
of the products of short distance coefficients and NRQCD long distance matrix elements. 
The short distance coefficients can be expanded as power series of as at the energy scale 
of the heavy quark masses, thus they are calculable perturbatively. We may attract them 
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by matching the results obtained with QCD full theory calculations at the threshold to the 
ones obtained by NRQCD effective theory. Generally, the long distance matrix elements 
can be estimated only by means of the velocity power counting rule of the NRQCD effective 
theory. 

The effective theory NRQCD considers a quarkonium state as an expansion of various 
Fock states. The annihilation may be via various Fock states and will depend on various 
long distance matrix elements. According to NRQCD, the annihilation Be light 
hadrons can be carried out via leading Fock states |c6i(^S'o) >, \cb^(^Si)g > and \d>s(}Pi)g > 
etc. With a naive order estimate, for the leading order the concerned annihilation may be 
the processes via the leading Fock state in velocity v but an order higher in QCD, i.e., 
(c6)i(^5'o) — » + gg, and those via a higher Fock state but an order lower in QCD 
comparatively, such as {cb)s(^Si) l^vi + g and {cb)^{}Pi) l^vi + g etc, because of 
the compensation for the Fock states and QCD couplings. Therefore, in the decay B^. — >■ 
light hadrons we should consider the various short distance processes. Potentially the 
annihilation can be used to probe the contributions from various components of the higher 
Fock states, and it is why we investigate the annihilation quantitatively in this paper in 
the framework of NRQCD factorization. Namely we estimate the decay rates accordingly 
by taking the potential model value for the color-singlet matrix element which relates to 
the wave function (derivative of the wave function) at origin directly and by means of the 
velocity scaling rule of NRQCD for the color-octet ones. 

In conventional potential model, a heavy quarkonium just is a bound state of the heavy 
quark and antiquark in color singlet 0]. Thus in the framework of potential model, when 
the annihilation of the meson B^. is estimated, only the color-singlet component is taken into 
account. As a result, only color-singlet matrix elements in NRQCD may be related to the 
wave functions obtained in potential model. 

In Ref. 0, some typical long distance contributions to the leptonic radiative decay have 
been taken into account, and results show that such contributions are negligible in compar- 
ison with the short distance contributions. In fact in the annihilation concerned here, there 
are also similar long-distance contributions. Therefore in this paper we also estimate the 
long distance contributions. 

The paper is organized as follows: The first section contributes to the introduction. In 
section II, we briefly outline the formula for the annihilation within NRQCD framework first. 
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then in subsection A. we present the method to calculate the short distance coefficients for 
color-singlet one and final result, which corresponds to the process (c6)i — > l^vigg\ and in 
subsection 5., we calculate the short distance coefficients for color-octet ones, which relates 
to (c6)8 l^^iQ of higher Fock states (color octet ones). In section 111, we estimate the 
typical long distance contributions of the color singlet leptonic decays using the so-called 
'generalized instantaneous approximation' |^. In section IV, we present the numerical 
results, a brief summary and discussions. 

II. NRQCD FACTORIZATION ANALYSIS 

In the non-relativistic heavy quark limit, there are several distinct energy scales in heavy 
quark and antiquark bound state system. The heavy quark masses ttIc and mb relate to the 
largest energy scale, while the off-shell 3-momentum mv and the typical binding energy mv"^ 
are the small ones, where v ~ 0{a,.,) is the relative velocity between the heavy quark and the 
antiquark. NRQCD effective theory is established by integrating out the energy scale effects 
at heavy quark mass. In the effective theory, the heavy quark and antiquark are described 
by non-relativistic two component fields. A physical state of meson can be decomposed 
into a sum of a set of Fock states: 

\B,) = 0{v^)\ch{'So)) + 0{v')\cU^Pr)g) + 0{v^)\cU'S{)g) + ... . 

According to the velocity scaling rule as indicated in the above expension, the probability 
of each Fock state scales as a definite power of v. Namely the leading Fock state of the is 
|c6i(^S'o)), whose probability is of order 0(f°). The next leading Fock states are \cb^{^Pi)g) 
and \ch^{^Si)g)^ whose probability is of order f^. In fact, 'certain order' in v comes from 
the integration for phase space, but here we restrict ourselves to consider the wave function 
expansion itself so we do not take into account those from relevant phase space integration. 
Throughout this paper we denote a state of (c6) pair with a subscript for color: 1 for color 
singlet and 8 for color-octet, but for its angular momentum quantum number we put in the 
parentheses accordingly. 

The annihilation of the c and the h quarks happens in short-distance comparing with 
perturbative QCD (pQCD) energy scale, because of c and h being heavy. The characteristic 
space-time of the annihilation is of order 0{l/m^ or 0(l/m{,). The size of the B^ meson is 
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of order 0{:^), where /i — nib "t^c/ ("^6 + "^c) is the reduce mass of the b and c system. They 
are distinctly separated before annihilating. Thus NRQCD factorization formula can also 
be used to analyze the decay of the Be meson. We adopt the formulism of NRQCD to do 
calculations in this paper without strong modification. In the annihilation Be light 
hadrons . the b, c quarks annihilate by weak and strong interactions both, that, according 
to pQCD, is of the short-distance effects, because it happens in such a short space-time at 
an order <^ i/^qcD- According to the NRQCD factorization formula, the decay width can 
be factored into 

r = ^^Cn{Be\dnW , (1) 

where C„(n = 1, 2, • • •) are the short distance coefficients, and {B^ On\Bc){n = 1, 2, • • •) are 
NRQCD long distance matrix elements. The short distance coefficients, being expanded in 
power of as, can be calculated with pQCD. While the long distance matrix elements are 
the 'average' values of the operators 0„, n = 1, 2, ■ ■ ■ which, being local gauge-invariant and 
relevant to the annihilation, consists of non-relativistic 4-quark operators. They measure the 
inclusive probability of finding a pair of cb being at the same point and with suitable color 
and angular-momentums being specified. By the the velocity scaling rule the possibility 
in the meson Be in order of v. The (be) pair in the meson Be, not only in color-singlet 
but also in color-octet, may be annihilated at short-distance. Correspondingly, the matrix 
element is regarded as a color-single one or a color-octet one. The lowest order process for 
the color- singlet, i.e. the short- distance annihilation of the meson Be to leptons and light- 
hadrons is (6c)i —>■ lui -\- gg, while those for the color-octet ones are {bc)^{^Pi) g and 

(6c) 8 (^5*1) — > Ivi -\- g. Hence the factorization reads as the following form: 

C,{'S^)^^{BeWcX, Ai^c\Be) 



2M 

+ CsCS,) {Be\4^'T'^XbxWT^HBe) 

+ CsCPi) ^ {Be\4Hl ^)T''^c\Be) ] , (2) 

where Ci(^S'o), Csi^Si) and Cs(^Pi) are dimensionless short distance coefficients and they 
are of order GpOi^, GpOLsi and GpOLs (the lowest order in weak interaction with suitable 
strong couplings) respectively. The matrix elements of the three term in the bracket scales 
as v^, and v^, respectively. Note that here we adopt the normalization for the meson 
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states as {Bc{P')\Bc(P)) = 2i?p(27r)^(5^(P' — P) and since we focus only those annihilations, 
in addition to the leptons, also to light hadrons, i.e. the pure leptonic annihilation of the 
meson is not included, hence the contributions from pure leptonic annihilation should 
not be included in the coefficient Ci(^Sq) appearing in Eq.(0). 

Now let us to compute these short-distance coefficients. Namely we calculate the processes 
with pQCD precisely, and by matching the results with the NRQCD factorized formula, 
finally we will obtain the coefficients. 

The matrix element for the decay Bc{P) J'^{k4) + i^iik^) + X can generally be written 

as 

M = ^ {lu,\m {X\J,\B^) , (3) 

where Vcb is the CKM matrix element, and j'^ and are the weak currents for lu leptons 
and cb quarks. 

The decay width can then be expressed as 

Gj^V], r d% d% 



J LEa 



where 



4M(27r)6 ] 2^4 Is 

t,.{P,k) = lmj d'xe''^B,{P)\j;{x)J:{0)\B,{P)) , (5) 

are the leptonic tensor and the hadronic one, respectively. The cb quark currents in the 
hadronic tensor precisely are = ipc'y^il — l^j'ipb and = ^/'b7i.(l — 75)V'c- 
The leptonic tensor L^'^ can easily be calculated and reads: 

L^- = S{ktkl + e,kl - g'^-'ik, ■ k,) - le^^'^'^^k^^k.^), (6) 

where £^^^'^l^ is the total antisymmetric tensor. 

According to the factorization, the hadronic tensor contains both short- distance coef- 
ficients and long-distance matrix elements for suitable operators. Based on NRQCD, the 
hadronic tensor can be factored as: 



V(P,fc) = \mj d''xe'^%B,{P)\J];:{x)j:mB,{P)) 



d^,^CSo;P,k)^{BMxb xl^c\B,) 
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+ d^CPu P, k) {BMH^ D)T\, xl{-^\ D)T^i^cW ] , (7) 

here di'^ Sq] P, k) , dg"^ {^Si; P, k) and dg'^ Pi, P, k) are the factors of the short- distance 
coefficients which we need to compute out precisely. Comparing Eq.(^) with Eq.(0), they 
are related to the coefficients Ci(}So),Cs{^Si) and Cs(}Pi) in Eq.(0) precisely as follows: 

Glvl rd%d^k 

Cs{^Pi) = 
with k = ki + k;j 



These coefficients can be evaluated by using the threshold expansion method[|T5|. In the 
method, on-shell (c6) pair near the threshold with specified quantum number are annihilated. 
The amplitudes for the processes at the quark level are calculated in the fundamental pQCD 
and expanded as the relative velocity v in the framework of factorization formula. In the 
meantime, the real process of the B^ annihilation may be calculated in the framework of 
the effective theory NRQCD and 'suffers' a similar factorization, in which the short distance 
coefficients are the same. By matching these two ways of calculations, the short distance 
coefficients can be 'read off'. 

More precisely the coefficients C7i(^S'o) and dl^i^So] P, k) are determined by matching the 
process (c6)i(i5o) ^ l+vi + gg, while Cs(^Si) and Cl^'^i; P, k)- Cs^Pi) and rffl^A; P, k) 
are determined by matching the processes (bc)s{^Si) l^ui + g; (hc)g(^Pi) l^ui + g 
with the calculations on the annihilation of the meson B^. with the effective theory NRQCD 
respectively. In the following two subsections, we show how to determine these coefficients 
by matching relevant processes. 



A. Short distance coefficient for color-singlet 

In this subsection, we use the threshold expansion method to determine the short distance 
coefficients Ci("'^S'o) and di'^(}So; P, k) by matching the process 0(^1)6(^2) ik4)ui{k^) + 
g{ki)g{k2), where (c6) is in a color-singlet state. At the lowest order, there are only six 
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Feynman diagrams contributing to the matrix element. Three of them are shown in Fig.|l| 
and the other threes can be obtained by exchanging the gluons. The amphtude is given by: 



where A^^ is defined by 



M = ^{lui\j,\0) , (11) 



A, = {gg\J,\{bc),CSo)) = E (12) 

j=l,2,-,6 



The six terms in the amphtude correspond to six Feynman diagrams, respectively. The first 
three, which correspond to the diagrams in Fig.|^, are given by: 

2s:ab _ 1 1 

A,,, ^ 6fe) -,,(1 - , (15) 

where b{p2), c{pi) are the Dirac four-component spinors of the antiquark b and the quark c 
respectively. e"'{ki) and £^(^2) are the polarization vectors of the two gluons with color in- 
dices a, b respectively. The other three terms of the amplitude can be obtained by exchanging 
the gluon vertices. 

The momenta of b, c quarks are related to their total and relative momenta P and q as 
follows: 

Pi = + q , P2 = fJ'bP - q (16) 

here Ur = — and Uf. = — . In the case of the lowest order calculations for the 
5— wave, that is we are considering, the relative momentum q can be set to null, i.e. q = 0, 
hence the momenta can be further simplified as 

Pi = fXcP , P2 = f^bP ■ 

Then the hadronic tensor t'^'^[P, k) for the annihilation [(c6) 1(^5*0)] ^ + i^i + g + g can 
be expressed as 

r'^(P^k) = ^ f^^A^^^A'^ 5\P-k-k,-k2) 



(27r)2 7 2ei 2e2 

d^,'CSo;P,k) {{cb),CSoMxi>xlMcb)iCSo)) + -- - , (17) 
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where "• • •" denotes those terms, which, as the ones corresponding to the pure leptonic decays 
and those corresponding to the (cb) with other quantum numbers, we are not interested in 
this paper. 

The integration over the gluon phase space can be preformed. The details of the calcu- 
lations are given in the Appendix A, but we outline the basic steps here. By introducing 
three Lorentz invariant variables: z = y = and Xi = ^^p^, the hadronic tensor 

t'^'^{P,k) can be expressed in terms of them. The short distance coefficients d^^ (}Sq] P,k) 
can then be obtained by matching the Eq.(|l^) and Eq.(|^). 

The most general form of the tensor di^i^So] P, k) can be expressed as: 

d'l^CSo;P,k) = Ag"" + BP^P" + CP''k'' + Dk^'P" + Ek^'k" . (18) 

However, when we contract it with the leptonic tensor Eq.(^), the contribution from the last 
three terms of Eq.(|l^) vanish due to the fact that we ignore the lepton masses totally, so 
the leptonic weak current (even including the axial component) are conserved. Namely we 
have k^L^'^lk^, k^) = 0. Thus only the first two terms are effective. The detail calculations 
and the expressions of the coefficients A and B are put in Appendix A. 

Since we are interested in the energy spectrum of the charged lepton which is measurable, 
so we also give the coefficients before carrying out the integration over the leptonic spectrum 
xi, and those before the integration over the z,y and xi: 

(~<2 T/2 yM'4 ^ 

CiCSo) = J dydzdxi{dt\^S,- P, k)L,,)^xi)^{l + y-z- xi)^{Y), (19) 

where the new variable 

^ {1 + y - z - xifxf _ _ {1 + y - z- xi)xi 2 
~ 4 2 ' 

is introduced. By setting Y = and drawing the Dalitz diagram for the process, we may 
determine the integration area for y, z and Ei (or xi). The result about the integration area 
is also put in Appendix. Since the integration is quite complicated so we carry out it (them) 
numerically. 

If one carry out the integration in the order: to integrate y and z first, then Xi, then, 
before doing the last integration for xi, the measurable energy spectrum of the charged 
lepton ^ can be obtained. Thus we will do the integration in this order and discuss the 
obtained energy spectrum in Section IV. 
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In order to check the numerical results, we also take a different order to do the integration. 
Namely we try to make the integration of the variable of the leptons y at the last, but try 
to integrate the other variables of the leptons and those of gluons. First of all, we need to 
change the lepton tensor from L^'^ to N^": 

TV*- = —{ki'k'' - k^g^"") . (20) 

To complete the phase space integration in this order, and having all the other variables, 
such as angles etc, of the phase space integrated out, we reach to the step that only three 
independent integration variables: y, the sum of the gluons' energies Es and the difference of 
the gluons' energies need to be integrated out. Then to carry out the integration further, 
we need to determine the integration area for these three variables by the Dalitz diagram, 
and the obtained result about the integration area is < |/ < 1, —£^°-^ < £d < ^max 
_^min <Ss< e^""^, where 

sT' = ^(M - SsY - M^y, 
So the coefficient Ci(^S'o) is written: 

/ max max 
dy r r fiies,ed,y)desded, (22) 
J^,ni,i J_^max 

where the integrand fi{es,ed,y) is obtained straightforwardly in this order step by step as 
described above. Since the results for each step are very tedious so we do not present 
here, but we have done the computations carefully and really have the check for the phase 
space integration numerically. Indeed we find that the final results for the annihilation 
rate obtained by numerical integration in these two orders are the same, so we are sure 
that the numerical results very well. In fact, as a semi- finished result in this integration 
order, the 'spectrum' V — be obtained too, but it is not easy to measure 

experimentally so we will not present the curves of the spectrum here. 



B. Short distance coefficients for color-octet 

In this subsection, we use the threshold expansion method to determine the short dis- 
tance coefficients Cs{^Si), dg^{^Si; P, k) and C^{}Pi), d'^^{^Pi] P, k) by matching the process 
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c{pi)b{p2) l'^ {ki)vi{k^) + g{ki) where cb is in a color-octet state obviously with the one 
computed by the effective theory NRQCD. At the lowest order, there are only two Feyn- 
man diagrams contributing to the amplitude of the process. They are shown in Fig.^. The 
amplitude is given by: 



M"' = ^(/^,|/|0X, (23) 



where A'^' is defined by 



K = Wipe),) = (24) 

1=1,2 

where a = 1, ... 8 are the color indices. The two terms of the amplitude, corresponding to 
the two Feynman diagrams shown in Fig.|^, are given by: 

<i = IttIT'^T^ b{p,) 7,(1 - 75) ^^ _ ^1 _ ^f {ki) c{p,) , (25) 

<2 = \Tr[T'^T']g,Kp2)f{h)-, j 7,(1 - 75) c(pi) . (26) 

In the cases of the S"— wave and P— wave cb that we are considering and to the lowest order 
approximation, we expand the expression in power of q, the relative momentum, and keep 
the terms only up-to linear ones of q (because we are doing the 'leading order calculations 
only). 

Then the hadronic tensor t'^'^{P,k) for the annihilation {cb)i{^So) ^^^W can be ex- 
pressed as 



+ d^'Pu P. k) ^ {{cblP,m-t]^ d)T^Xb xIhI D)T^^c\{cblP,) ] + ■ ■{27) 

where ■ ■ ■ denotes the terms corresponding to the cb being in the other states, which we are 
not interested in here. 

The integration over the gluon phase space can be performed easily. Thus with Eq.(^ 
and Eq.(0), the short distance coefficients can be computed easily. 



C^i^S^) = / dxi^ 2 2 - 4 1 - ^i) log 1 - xi)) 28 
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for the component ((c6)8(^5'i)), and 
C8( Pi) = / dxi ^ —\xi 



M^{ - ?,mbml{xi - 1) + ml{xi - if + ml{xi - 1) 



127r^mlm^{l — xi 

- mlnicixi - l){8xi - 7) + 2m^m^[3 + 2{xi - 2)xi]} 

+ 2Mmbmc{ - ml{xi - l){2xi - 7) + m^mlixi - 1)(1 + 2xi) 
+ m^[(7 - 2x])xi - 5] + mb^mc[l + Xi{^Xi - 11)]} 
+ 2mlml{ml[ll + {xi - Q)xi] + mlKxi - 2)xi - 1] + 2mbmc{x^ - 3)} 

- (1 - xi) ln(l - xi) M'^{mlmc{7 - 6xi) + ml{xi - 1) + ml{xi - 1) 
+ 2mlml{2xi - 3) + mbml{2xi - 3)} 

- 2M'mbmc{5ml{xi — 1) + nibml^xi — 1) + mlrricil — 7xi) + ml{5xi — 7)} 

+ 2mlml{ml{xi + 1) - 2mbmc{xi - 3) + m^(5x/ - 11)} | (29) 

for the component {{d>)s(}Pi)). 

III. THE LONG DISTANCE EFFECTS 

Besides those between b and c quark pair in the initial state (the interactions make the 
two quarks into a bound state:the meson Be), there are also multi soft-gluon interactions 
between the c and c pair when the 6-quark has decayed into c and the pair of c and c has 
not annihilated yet. Since the interactions owes to multi soft-gluon exchange, so people 
generally attribute them as to long distance effects. In this section, we consider this kind 
of effects in the annihilation for the color singlet component Be — ^ l^vigg only, namely the 
possible ones, where the interactions make the cc pair to form a bound state rjc with suitable 
quantum numbers, because they may be sizable based on naive order counting but for the 
concerned octet components the effects cannot be very great. 

As a rough estimate, it is enough to consider the dominant effects just to take into account 
the possible intermediate meson state rj^. Now what we are considering is described in Fig.|^. 
The amplitude corresponding to Fig.^ is 



2V6 ' ' ' " /^vo,,v .v,,^ ^^^^ 

d^r ^ , ^ 1 , „ . _ 1 



Tr[xQ{r)-i^ ^ 7„] + Tr[xQ(r)7n ^ 7„ 



(27r)4 V ^"^^ "'^\Q + i-h-^c' ' '"'hQ + f-h-^c 
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■ / (0i^^[^^(^)rM^Q(^')(/^c/' + i + mj] , (30) 

where V ^ = 7^(1 —75), £^(^1) e'^{k2) are the polarization vectors of the two real gluons, 
q = yUfePi— /icP2 and q' = q+ ^[{fic — fit) P + k] and r = ^{pi—p2 + k) are the relative momenta 
between the two constitute quarks of the Be meson and the bound state rjc correspondingly. 
Note that the equations 

P = P1+P2 , k = ki + k^ , Q = ki + k2 



have been used in Eq . (|30D . uJq = yQ'^ + (M')^ (M' is the mass of rjc). Note that for 
the time-component, usually off mass shell we have Qq 7^ u^, but in the present case, the 
charged lepton energy may be not very high, i.e., Ei < ~ ^ , so the intermediate bound 
state ?7c may reach to its mass-shell, that is Qo = and leads to a singularity in Eq. (pO|) . 
This difficulty happens is due to the fact that we ignore the width of 7]^. Thus here the 
width of r]c should be considered i.e. to replace ie with ^-^^^ in the relevant propagator, 
where ie is the 'tiny' imaginary quantity but V^^ is total width of the bound state rjc of cc. 
Then under the non-relativstic approximation, we have 

1 1 



Since we consider rjc as the only one intermediate bound state as indicated by Eq. (pO|) , 
so only the 'weak current matrix element' corresponding to Fig.|^ and the amplitude for rjc 
annihilation are needed to be computed: 

{r,^T,\P) = z J J^J^r[xBAq)^,X,imi + ^c)] , (31) 



and 



+ Tr[x,.(r)7.^^3^7j) • (32) 

In fact it is easy to compute the second matrix elements, while to compute the first one, 
that corresponds to the Feynman diagram Fig.|^, we need to pay more care to take into 
account the effects of recoil for the intermediate state. We adopt the approach, i.e. the 
so-called generalized instantaneous approximation, which was proposed firstly in Ref. [Q, 
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to deal with the recoil effects in the first current matrix element. The main points of 
the approach are to 'extend' the potential model, which is based on Schrodinger equation, 
into the one on Bethe-Salpeter (B.S.) equation for the non-relativistic binding systems, 
and then, according to Mandelstam method []TB|, to formulate the current matrix element. 



so the current matrix element is set on a full relativistic formulation, finally by making 
the 'generalized instantaneous approximation' on the whole relativistic matrix element, i.e., 
integrating out the 'time' component of the relative momentum in the formulation of the 
matrix element by a contour integration, the matrix element turns back to be formulated by 
means of the Schrodinger wave functions sandwiched by proper operators. The Schrodinger 
wave functions are just those of the original potential model for each system and have 
direct relation to the B.S. wave functions through the original potential model for each 
system and have direct relation to the B.S. wave functions through the original instantaneous 
approximation proposed by Salpeter. Since the approach at a middle stage has a fully 
relativistic formulation for the weak current matrix (the Mandelstam formulation of the 
matrix element), so the final formulation surely takes the recoil effects into account properly. 
Since one may find the details of the approximation in several references[|, ^ ||, and so we 
do not repeat it here. 

By applying the result of the generalized instantaneous approximation on the current 
matrix element to computing the matrix elements, one may obtain the final formula straight- 
forwardly and estimate the long distance effects. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 



We take the parameters as in Refs.[|T], |, 0: a, = 0.24, Tb, = 2.7Ups~\ ip{0) = 

3 5 

O.SSOGeV^a, ip'i^) = 0.250GeV2 ^ in numerical calculations. For comparison, let us list here 
the branching ratios of the pure leptonic decays for Be meson 

Br{Be eue) = 1-89 x 10"^ , 5r(fi, ^ fiu^) = 7.57 x 10"^ , 



Br{Be TUr) = 1.95 x 10"^ . 

To obtain more reliable values, we should be careful to take the quarks' masses. According 
to the discussions in Ref. |Q, here we take the effective masses of c and b quarks to be 
m^'^'^ = 1.5Gel^, m^-^-^ = A.QGeV . For the mass of B^, meson, we take pole masses of the 
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quarks m^°'^ = 1.88GeV, m^°^'^ = 5.02 [0, 0, [T^ first, and then being consistent with that 
of potential model, the value of Be mass M = 6.352GeV is taken. 

For the short distance contributions of the process —>■ gglv-, we find that when Ei < 

2 2 

— where Ei is the energy of the charged lepton, being very different from that of the 



2mi, 

one photon radiative correction 0, the c quark from the decay of h may reach to mass-shell. 
Thus to obtain meaningful result, here we should also keep the 'width' of c quark in its 
propagator, namely we should make the following replace: 

i + nic i + rric 



q'^ — ml + it q'^ — ml + i{4:mcTc) ' 

where is the total width of an on-shell c-quark and from D decays its value should be 
Fc ~ 1.229ps-i §. 

For the annihilation of color singlet component. Be —>■ gglv-, since the color singlet matrix 
element {Bc\ip\xb xIV^cl-Sc) may be related to the wave function at origin squared ['0(0) p 
and if only the short distance contributions are taken into account, the branching ratio can 
be computed out quite precisely: 

Br'^°^\Be luigg) = 2.71 x 10"^ . (33) 

Whereas if only the long distance contributions are concerned as in Sec. Ill, the branching 
ratio may be computed out too: 

5r'°"3(Bc ^ luigg) = 4.45 x 10"^ . (34) 

From the above values one may see that the long distance effects are quite great here. Note 
that there is slight overlapping for the short distance contributions and the long distance 
contributions and these two kinds of contributions should have interference, but here we 
ignore all of them. 

For the annihilation of color octet components, B^ li^g, there is no reliable way to 
calculate the color octet matrix elements {Bc\ipla'^T"'Xb Xb'^^T""ipc\Bc) and {Bc\ipl{—i^ D 
)T°'Xh Xfel""^! D)T""ipc\Bc) which are necessary when computing the annihilation. 

In order to have order estimate roughly, we, based on the velocity scale rules of NRQCD, 
try to assume them to be smaller by certain order 0{v) than the S-wave wave functions 
at the origin |?/^(0)p for color singlet, and the derivative of the P-wave function at origin 
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|'0'(O)p, which can be computed by potential model. Namely based on the velocity scale 
rule, we assume 

and 

{bMhI D)T'^Xb xI(-4 B)r>,|s,) ^ Al ■ {b,\4hI D)xb xl^l omB,) 

with As', Ap being 0{v) constants. So with the assumption we may evaluate the color octet 
contributions accordingly. 

In order to explore the characteristics of the color octet components in the decay for the 
meson Be and to have a comparison with those of the color singlet, we try two possible 
choices for the S-wave color octet wave function and the derivative of the P-wave color octet 
wave function at the origin: CASE A: A5 ~ Ap ~ 0.1 and CASE B: A5 ~ Ap ~ 0.3, 
because we think that according to the velocity scale rule the range from CASE A to 
CASE B is reasonable and it is helpful to see the possibility if the color octet components 
of the meson can be detectable or not experimentally. The branching ratio of the color 
octet component annihilations for CASE A: 

Br{Bc{ch{^Si)) luig) = 1.73 x 10"^ , Br{Bc{ch{^Pi)) Ivig) = 2.24 x 10"^ (35) 

and for CASE B: 

BriB.iSsCS,)) ^ luig) = 1.55 x lO'^^ , Br{B,{cbsCPi)) ^ luig) = 2.02 x 10"^ . (36) 

From the values above, we see that the helicity suppression in the annihilation processes 
of Be is released. In total, the color singlet annihilation to light hadrons is bigger than 
those of color octet. Furthermore, when there is one more gluon bremsstrahlung than the 
pure leptonic decay, there are more freedoms in energy and quantum number carried by 
the gluon(s), that the intermediate, c quark, even a relevant bound state such as rjc etc 
(compounded by the produced c-quark and the 'original' c-quark in the meson Be) may 
become on shell, thus the width of the concerned annihilations may be so great as that of 
the semi-leptonic decays b civ or B^ riJu correspondingly. 

In Fig. ^, we show the lepton energy spectrum of the annihilation to leptons and light 
hadrons. The dashed line dictates that of the spectrum for the color singlet components, 
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where the short distance and long distance contributions are combined. The dotted and sohd 
hnes, which correspond to the contributions from the color octet components (c6)8('^>S'i)) and 
(c6)8(^-Pi)) in Be meson respectively, dictate the spectra for the color octet components. It 
is interesting to point out from Fig. |^ that because the meson Be the wave function of the 
color octet components are suppressed, the concerned annihilations due to the color octet 
components (one-gluon bremsstrahlung) even though are smaller comparatively in the most 
region than that due to the color singlet component (two-gluon bremsstrahlung), it may 
become greater in certain region of the spectrum, so that it is possible to see the color octet 
contributions experimentally though studying the changed lepton energy spectrum of the 
inclusive decay Bf. ^ l^vr ■ ■ carefully, especially, around the end point of spectrum. 

In order to describe the difference quantitatively in the spectrum of the charged lepton 
near the end point, where the color octet contributions may become dominant over those of 
the color singlet ones, we introduce the ratios of the integrated partial decay widths of Be 
meson for the color-singlet and the color-octet, 

' r(c6i(i5o),xr*) + r(c68(3^i),a;r*) 
for the S-wave color-octet component |c68(^'S'i)), and 

r(c&i(igo),xrO 

T{ch{}So),xr') + T{ch{^P{),xr') 

for P-wave color-octet state |c68(^-Pi)), which depend on the cut of the lepton energy x^*. 
Here 

ax 
dx 

and 

Jx'P* dx 

with 5 = and a giving tiny gluon mass nig = 0.2GeV. We evaluate them and put result 
in the Table I. 

From Fig. ^ and Table I, one may see clearly that there is a possibility to verify experi- 
mentally if the color octet components in the meson Be play some roles in the Be annihilation 
Be —>■ l^ui+hadrons. It is known that NRQCD is a very absorbing theory, whereas it needs 
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to verify widely thus we thank that our estimates are very prehminary but may addresses 
experimentahsts' attention on this subject. It is worth further to study the possibihty in the 
hadronic coUision environments more quantitatively and the first is Monte Carlo simulation 
and then carrying on experimental analysis for the verification. 
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APPENDIX: INTEGRATION FORMULA OF THE COLOR SINGLET Be MESON 

A number of kinematic variables will appear repeatedly in the discussion. For clarity all 
of them will be collected first. P, ki, /c2, k^, k^ denote the four- momenta of the various 
particles, Q = ki + k2, and k = k^ + k^ characterize the gluon-gluon system and the virtual 
W respectively. 



The scaled masses and lepton energies 
vary in the region 



y=T7^.^ ^ = 77^' ^' = irr' = ^ (37) 



< < 1 , 
< y <xi, 

0< z < w = (l-a;0(l--), (38) 

Xi 

where the leptons' mass is ignored. Frequently used kinematical variables which characterize 
the gluon-gluon system are 

jl-y + z) 
-n-o — ) 



Rs - 2 ' 

Yp = In (^(^i±M^ . (39) 

All of the variables are chosen in such a scaled form on the Sc-meson's mass so as to make 
no explicit mass dimension left for convenience in the following defined coefficients q. 
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To calculate the annihilation, one observes that the squared matrix element can be 'fac- 
torized' into the leptonic tensor L^y{ki, k^) and the hadronic tensor T^i,{P, ki, /C2) = A*^ ■ A^, 
(A^ is the amplitude). To perform integration over the gluon-gluon space, the resultant 
integration t^^, depends on P and Q = ki + k2 only: 



V(P,Q) = I dR2{Q;kuk2)T^AP.ki,k2) 



= rfi,M.('5o; P, k) ^ {BMxb xl^cm (40) 
where the phase space dR2 is defined by: 

dR2(Q; fci, k2) = (27r)45(Q - k^ - k2)-^ /'I' , (41) 

and the tensor di^f^^{^So; P, k) has the general structure: 

di,^XSo\ P, k) = Ag^, + BP^P, + CQ^Q, + DP^Q, + EQ^P, . (42) 

For the massless lepton pair using the lepton pair current conservation relation, k^L^^ = 
k^L^^ = where k = P — Q is the total momentum of the lepton pair. With these equations 
the above equation Eq.(^2D can be further simplified as 

di,^,CSo; P, k) = Ag^, + (B + C + D + E)P^P, = Ag^, + B'P^P,. (43) 

Instead of integrating the tensor T^^^, it is sufficient to integrate the following scalar 
projections 

ci = / dR2T.,P^g^'' , 

C2 = 77^^ / dR2T^,P^P'', (44) 



and the coefficients A, B', ■ ■ ■ can be expressed in terms of the scalar functions q, i = 1, 2 as 
follows: 

A='-^, B' = ^^^. (45) 
3P2 ' 3p4 y ^ 

The numerators of Cj are given by polynomials in (Pki) and {Pk2) with coefficients 

depending on y, (QP) = M'^Rq and (Qki) = {Qk2) = In fact, all of the phase space 

integrations may be attributed to the following integrations: 

= j dR2{Q;k,,k2){Pk^nPk2r (-2 < m,n < 0). (46) 
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Because of the Lorentz invariance, we evaluate the integral in the Q — ki + k2 rest system 
and the results show 



-^0,0 



1 

1-2,0 — Io,-2 — 



ATrM<^^RsRl zRq' 

Yp 



47rM4i?oi?3' 

1 Y 2 

^-^■-^ ^ I^^R^, + JWj- ^^^^ 
The explicit form of Cj are shown in the following, where a common factor ^o;^ is con- 
tracted out for convenience and the terms that can be obtained by interchanging rrib and rric 
is not shown explicitly, that is, the actual value of each Cj equals (q + Ci(ri <-> r2)), where 
ri - ^ and rs - ^ 
1 fi 

ci = ,2 , 2 2 p2 (-8^^2(3/2^ + 4(/2 - A)r,^)i^g + f^^yz + 4it:g) 

-4/irir2(/i/2^ + 2/2(/i - z)R, + (2/1/2 + ^ - 4/2^ + 4/2r2)i?^) + 
8rl{f2z{-fi + z)Rl + f2T2Rl{'iz + (4/1 - 2z)Ro) - 4f2r2^Ro{2f, + Ro + R^) + 
rl{2flf2 + i?^(4/2 - 4/1/2 + 5/iz - 2/ii?o + 4/2i?^)))) + 

-{lQf2zr\Rl + flf2z{y -2Rl) + 



/I r J 4d2 

nf2rlrlRlR, 



1-3 

4/iri(-(/i/2;2r2) - 2/2(/i - z)r2Ro + r2(2/i/2 + z{z - I - 2/2) 

+4/2(3 - 2r2)r2)i?^ + 2{f2Z - 2zt2 + V2r''M + 4r2i?o) + 8r^2i?o(3/2-2 + 

4r2(/ir2 + hRo)) + 8r?(-2/2zi?2 + 

4/2r2'i?^ + /2(-2/i + z)Rl - 4/2r2'i?o(2/i + Rl) + 

/2r2i?^(^ - 2i?o(2 - 2/1 + z + 2i?^)) + ri(2/iV2 + 

i?o(-4/i/2 + 4/iz + 3/2^ - 4(/i - /2)i?o(l + i?o))))), (48) 

C2 = f2fj^2^2ji2 ifiM^' + 8rir2(-^ + 4rir2)) - 

4/i/2rii?o(r2(4/i - ^i?o) + ri(^i?o + 16r^(-2r2 + i?o))) 
+2Rl{flf2Z + 2ri(/i/2^ - 2/ir2(^ + 3/2Z + 4/2r2) - 
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Ary.ih + 2(-/i + /2)r2 - 2hRo) + n{hz' + 
2rl{hhz + 2/2(1/ + z) + 2hr^{l + 2r2 - 2i?o) - 4/ii?o))))) + 
f2f 2 Lap (/i/2^^ + 8fj2rM-z + 4rir2) - 16/i/2rir2(/i - S^r^^ji^o + 
2{flh{l - 2z) + 2ri(/2z(/i + (-2 - 3/i + Ay)n) + h{Ahh + 
(-3 + /2 + |/)^)r2 - 4(2/1/2 + ri (8/1/2 - (3/1 + f2)z + 
2/i(-l + ri)ri))r2' + 8/ir^2' + lQhrir2''))Rl - 
16ri(-(/ir2(l + f2-y + 2hr2)) - rl{f2Z + r2(/2 - 
(/i - 2/2)r2)) + rir2(/2(-4/i + y + z) + 

r2(2/i + h + (/i + ^f2)r2)))Rl - 16riV2(2/2ri + (/i - 3/2)r2)i?^), (49) 

where /i = -1 + y + 2rii?o, /2 = -1 + y + 2r2i?o. 

After contraction of the hadronic tensor with Iv/ii/(^4, ^5), replacing Q hy P — k and 
substituting 

(Py.^^,(P,).Ml+l^. (50) 

= 1 +y - 2; - X;. (51) 

One is left with the task to integrate the function of x;, y and z, which can be done numer- 
ically. 
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TABLE I: The ratios of the integrated partial width Rs and Rp for CASE A {As = 0.1, Ap 
and CASE B {As = 0.3, Ap = 0.3) (The definition about Rs, Rp, As, Ap and x^* is in 





CASE A: {As = 0.1) 


CASE B: {As = 0.3) 


cut 


0.80 


0.85 


0.90 


0.95 


0.80 


0.85 


0.90 


0.95 


Rs 


0.257 


0.191 


0.122 


0.0-17 


0.037 


0.026 


0.015 


0.005 




CASE A: (Ap = 0.1) 


CASE B: (Ap = 0.3) 


cut 


0.90 


0.92 


0.94 


0.96 


0.80 


0.85 


0.90 


0.95 


Rp 


0.280 


0.212 


0.154 


0.089 


0.132 


0.081 


0.041 


0.014 
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FIG. 1: Three typical Feynman diagrams for the decay Be l^vgg- 



24 





FIG. 2: The diagrams for the annihilation width of Be I'^i^g, where Be is indicated in a color 
octet Fock state. The vertical curcul line in the diagrams is understood the according imaginary 
part being taken. 
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FIG. 3: The diagrams for long distance effects. 
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FIG. 4: Current matrix element for the long distance effects. 
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FIG. 5: The energy spectra of the charged lepton with different 'color-octet matrix elements'. The 
left figure is that for CASE A: A5 = 0.1 and Ap = 0.1. The right figure is that for CASE B: 
A5 = 0.3 and Ap = 0.3. The dashed line in the figures stands for the summed effects of short 
distance contributions and the long distance one with r/c as an intermediate bound state in the 
color singlet decay I'^vgg. The dotted and the solid lines for the color octet annihilations 

Be I'^J^g with cbsi^Si) and cbs(^Pi) components in Be meson correspondingly. 
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